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Abstract 

Cumulant expansion is used to derive accurate closed-form approx- 
imation for Monthly Sum Options in case of constant volatility model. 
Payoff of Monthly Sum Option is based on sum of N caped (and probably 
floored) returns. It is noticed, that 1/y/N can be used as a small parame- 
ter in Edgeworth expansion. First two leading terms of this expansion are 
calculated here. It is shown that the suggest closed-form approximation is 
in a good agreement with numerical results for typical mode parameters. 

1 Introduction 

Monthly Sum (MS) Options are embedded into popular fixed index annuities. 
In contrast to plain vanilla options, the closed-form formula for MS derivative 
is not known even in the simplest case of constant volatility and we need to use 
Monte-Carlo simulations to price them. To derive accurate approximation for 
MS options we can use some sort of expansion with small parameter. In the 
next section we will shown that MS derivative can be considered as a Call option 
on underlying with caped (and probably floored) returns. Caped distribution is 
very different from Gaussian one but, according to the central theorem, in the 
limit of large number of returns (N —> oo), the final distribution at time of expi- 
ration will be close to the normal one. So, in the leading 1/y/N approximation, 
we can use standard European Call formula to price MS options [1], [2]. 

To calculate 1/y/N corrections we need to choose convenient way to approx- 
imate non-normal distributions of caped underlyings. Various types of expan- 
sions in series of Hermit polynomials are used in science to approximate dis- 
tributions: namely Gram-Charlier, Gauss-Hcrmite and Edgeworth. But Gram- 
Charlier series has poor convergence properties [3]. Gauss-Hermit expansion 
is better but it has no expansion parameter and therefore it has no intrinsic 
measure of accuracy. The most prominent is Edgeworth expansion which has 
an expansion parameter l/\f~N [4]. This expansion is a true asymptotic series 
and has good convergence properties. 
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Here we consider Black-Scholes case of constant market volatility. It is 
shown, that for typical Monthly Sum option and market parameters, first two 
leading terms of this expansion give a good approximation for MS derivative. 
Note, that these first two terms of cumulant expansion are identical for all three 
types of cumulant expansion. But next corrections need to be grouped as terms 
of Edgeworth series to be a true asymptotic expansion with 1/y/N as a small 
parameter. Higher terms of cumulant expansion can improve accuracy of this 
approximation. 



2 Monthly Sum Derivative 



Monthly Sum option (MS) is defined as a derivative which payoff at the time of 
expiration t n = T is 



where 



(MS Payoff){cap) = Max ^ r m (cap), O^J 



r m (cap) = Min (S(t m )/S(t m -i) ~ 1, cap) 



(1) 



(2) 



are caped returns. 

In the case of small returns we can use logarithmic returns to approximate 
returns (2) 



r m (cap) ~ r m (cap) = Min Qn(S(t m )/S(t m -i)),ca$) 

where cap = ln(l + cap), and the derivative payoff is 

(MS Payoff)(cap) ~ (MSLN Payoff)(cap) 

~S(t N ) 



(3) 



Max 



N 



cxp ^2 r m (cap) - 1,0 



= Max 



S(t ) 



1,0 



where 



S(t N ) = S(T )e^ m =i r " 



(cap) 



(4) 



(5) 



It means that we can consider MS option as a call option on caped distribution 
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MS ~ MSLN = —Call*(X = S(t ),t N ) 



(6) 



where Call*(X) is a call payoff on the caped underlying 5(tjv) 

Fig.l demonstrates that, for typical market parameters, MSLN can be used 
as an accurate approximation for Monthly Sum derivative. 
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Figure 1: One Year MS and MSLN derivative values. Volatility=20%, Rate=3%, 
Dividend Yield=2%, Term=l year. 100,000 Monte-Carlo scenarios. 

3 Distribution Function and Cumulant Expan- 
sion 

To calculate call option value in eq.(6) and therefore to calculate MS derivative, 
we need to know distribution of caped underlying S(T) (5). 

To approximate this distribution at expiration time, it is convenient to repre- 
sent probability distribution function of logarithmic returns at time of expiration 
T in the form of infinite series as 

1 (*-^) 2 / fx - vT\ 

fx — uT\ fx — vT\ , . 

-*( wr) + H wr) + - ; (7) 

where x = \il(S(T)/S(0) is a logarithmic return, S(T) is a price of caped un- 
derlying at time T, ipn{x) is n-term of Edgeworth expansion, and H n (z) are 
Hermit polynomials: 

H n (z) = (-l)»e« /2 ^e- z / 2 ; H 3 (z) = z 3 - 3z. (8) 

Therefore the value of Monthly Sum derivative is 

MS ~ MSLN = M5 (0) + MS {1) + ... (9) 

where 

POO 

MS (n) = e -rT / ( e x _ y ^ n ^ dx ( 1Q ) 
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are terms of Edgeworth expansion for Monthly Sum option, r is interest rate. 

Calculating moments of distribution function ip we can determine parameters 
of this expansion: 



vT = «i; v 2 T = k 2 ; 6ei\/(v 2 T) 3 = n 3 ; . . . (11) 

where n n are cumulants. 

Using the fact that cumulant of sum of independent variables is equal to 
sum of cumulants, it is easy to calculate cumulants of longer term distributions 

as 

K n = Nb n ; (12) 

where N is number of time intervals (months) and i n is a cumulant of caped 
monthly distribution. 

Monthly cumulants can be expressed in terms of moments as 

H = Ii; i2=h-lh Lz=h-2>hh + 2ll; ... (13) 

where 

f cap x n dx (g-^) 2 f°° dx (x-^dt) 2 

I n = / e 2« 2 At + (cap) 11 / e ^a* (14) 

J-oc. V2ira 2 At J C a P V2ira 2 At 

are moments of caped distribution, At — t„ — t n —\ = 1/12 is a time interval, a 
is a market volatility, 

(j, = r - y - id 2 (15) 

is a risk neutral drift, y is a dividend yield. 

First two moments of caped distributions 1\ and I2 give us parameters of 
leading order approximation for distribution function ?/>(°) of caped underlying: 

vT = NIi ; v 2 T = N(I 2 - Ii)- (16) 

It means, that in the leading approximation, MS derivative value is an European 
ATM call option with volatility v and drift v. To set correct value for drift we 
need to adjust dividend yield y to its effective value 

y = r-v-\ 2 (17) 

Then, the price of Monthly Sum derivative in the leading order of cumulant 
expansion is 

MS^ = ^ Call(X + S(0),v, r, y), (18) 

On Fig. 2 we can observe that the accuracy of leading order approximation 
for One Year MS option with cap = 2.5% is not good for market volatility 
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Figure 2: One Year MS Derivative Value. Cap=2.5%, Rate=3%, Div Yield=2%, 
Term=l year. 100,000 scenarios were used in Monte-Carlo calculations (MC). 



a > 10%. To improve accuracy we need to take into account higher terms of 
cumulant expansion (7). 

The first correction to the leading term of cumulant expansion for probability 
distribution function (7) is 

1 (x— , T f (x — vT\ 
^ x ) = e 1 ^=L=e- L ^H 3 (^-^) ] (19) 



V2irv 2 T V vVT 

and according to eq.(13) 

£1 = N ^-^^ (20 ) 

As can we see on Fig. 2 this correction significantly improve accuracy of approx- 
imation. 

Formulas for moments of caped distribution and integral for (n = 1 

eq.(10)) are presented in Appendix. A. 

We also can consider MS derivative for caped returns with floor: 

r m (cap, floor) = Max (Min (S(t m )/S(t m -i) — 1, cap) , floor) (21) 

In this case, moments of this distributions can be calculated as 



/„ = (floorf 



fl0 ° r dx (*-Mt)2 

P 2ct 2 At -\- I — P 2cr z At 

J floor 

+ (cap) n / - e—^KT (22) 

J cap v / 2w 2 A< 
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First three moments of this distribution are calculated in Appendix. A eqs.(A.12- 
A.14). 



4 Conclusions 

In present paper we apply Cumulant Expansion to derive closed form approx- 
imation for Monthly Sum options. It is shown, that for typical MS derivative 
and market parameters this approximation works well. Calculation of higher 
terms of Edgeworth expansion can improve accuracy of derived approximation. 
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Appendices 



A Moments of caped distributions and correc- 
tion to MS price. 

Below, we present formulas for the first three moments of caped distribution 

/oo 

x n </>(x)dx (A.l) 
-oo 

where 

1 (x-nAt) 2 

4>(x) = 6 (cap - x) - p !» 2 4t + CS(x - cap). (A.2) 

V 27T(T 2 At 

C = r - dZ = 1 - CND(c); 

J cap V2w 2 At 

cap — /iAt 



CND(z) = f ^e-^ 2 = ^ (l + er/(s/V2) 



. 2tt 

e w = { i if,; (A-3) 

and <5(x) is Dirac <5- function. 

/i = -a-e"5 5 %/-i + /x AtCND(c) + C cap (A.4) 

V Z7T 

/ 2 = a 2 AtCND(£)--^a 2 Ate- £2 -2fiaAte~ £2 / 2 \ ^ + 

V27T V 27T 



(ji£tyCND(c) + C* (cap) 2 (A.5) 



+ 3^ At |a 2 Ai CND(c) - (cap - p At)e" £2/2 y j - 

30* Ai) 2 cr e" £2/2 + (p Ai) 3 CWL>(c~) + C (cap) 3 (A.6) 

V 27T 

First correction to the leading approximation (10) has integral which can be 
presented in the following form: 

_ v)e -^f + ( v 2 T) 3/2 e (.+ iv 2 )T CND ((^ + «) VT)(A.7) 
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To calculate MS derivative value with local floor, we need to calculate inte- 
grals in eq.(22) . This integrals can be presented in the following form 

I n = {floor) n CND(f) + (cap) n (l - CND(c)) + 

f cap x n dx /k 
+ / e 2^a* (A.8) 

J floor v27TCT 2 At 



where 



- = floor -/.At 
V2ira 2 At 



Using that 



ca P x n dx c-M')' (a 2 At)"/ 2 _ x 2 /2j ^ A 

e = ^ — / (x + /i)"e x /2 dx (A.10) 

/Zoor V27TCr 2 At V27T J/ 



where 



/i = -VAi (A.ll) 



we obtain the following formulas for the first three moments of the caped-floorcd 
distribution 

h = a (e~ f2 / 2 - e' £2/2 ) J ^ + fxAt(CND(c) - CND(f)) + 
V / V 2tt 

+ floor CND(f) + cap (1 - CND(c)) (A.12) 

2 At 

/2^r 



^A^CWD® - CND(f)) - ^ (ce" 2 */ 2 - /V /2 / 2 



. / /At 



2tt 

(MAt) 2 (C7VL>(£) - CND(f)) + 

(cap) 2 (l - CND(c)) + (floor) 2 'CND(f) (A.13) 



(cr 2 At) 3 



-((2 + ~c 2 )e-?/ 2 -(2 + f 2 )e-i 2 / 2 )sJ 
Z^aAt) 2 (cND(c) - CND(f) - (cxT^ >/2 - fe~^' 2 )^) 
3(^At) 2 cr (e 



e -c 2 /2 _ g-/ 2 / 2 



'At 
2^ + 

{nAtf{CND{c) - CND(f)) + 

(cap) 3 (l - CND{c)) + (floor) 3 CND(f) (A.14) 
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